1 



o 
o 



Y. R. Shimizu a , M. Matsuo b and K. Yoshida c 

a Department of Physics, Faculty of Sciences, Kyushu University, Fukuoka 812-8581, 
Japan 



Theoretical study of the decay-out spin of superdeformed bands in the 
Dy and Hg regions 

(N 

< 

o ; 

b Graduate School of Science and Technology, Niigata University, Niigata 950-2101, 
Japan 

>. 

c Institute of Natural Science, Nara University, Nara 631-8502, Japan 

qq ■ Decay of the superdeformed bands have been studied mainly concentrating upon the 
decay-out spin, which is sensitive to the tunneling probability between the super- and 
normal-deformed wells. Although the basic features are well understood by the calcula- 
tions, it is difficult to precisely reproduce the decay-out spins in some cases. Comparison 
of the systematic calculations with experimental data reveals that values of the calculated 
decay-out spins scatter more broadly around the average value in both the A ps 150 and 
190 regions, which reflects the variety of calculated tunneling probability in each band. 
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1. MODEL OF THE DECAY OF SUPERDEFORMED BANDS 

Decay of superdeformed (SD) rotational bands out to normal deformed (ND) states 
is one of the most interesting nuclear structure problems. It can be viewed as a shape- 
coexistence phenomena, and the mixing mechanism of two states having very different in- 
ternal structures can be studied as functions of angular momenta and excitation energies. 
Using the compound mixing model 0, we have investigated the decay-out phenomena 
in |f2[, where a consistent description was presented and the rapid decay-out was well 
understood. After the calculation of , more systematic and improved calculations have 
been performed || , in which all the relevant quantities to the decay-out have been calcu- 
lated without any adjustable parameters. We have found, however, it is rather difficult 
to reproduce the decay-out spin for individual SD bands (cf. Figs. [I], @). 

The basic idea of our framework |],[| is based on a simple two-well-mixing: The two 
groups of unperturbed states belonging to wells, I (e.g. ND) and II (e.g. SD), are separated 
by a potential barrier in some collective coordinate space (e.g. deformations). The most 
important is the (mean) coupling strength between the two, and is estimated as ||] 



J n-i 



(HuJi/27r)(huJii/27r)Tii-i, (regular-regular), (la) 
(Di/27r)(hu>ii/2ir)Tn-i, (chaos-regular), (lb) 
(Di/27r)(Dn/27r) Tn-i, (chaos-chaos), (lc) 



where ujji is the frequency of each well, Din = 1 / pi,n is the mean level distance, and 
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Figure 1. 7-ray intensities as functions of spin for nuclei in the A « 150 region. Solid, 
dashed, dotted and dash-dotted curves denote calculated ones for bands with different 
parity and signature sorted by excitation energy in the feeding region. Circle, triangle, 



square and cross denote experimental ones for yrast and excited bands taken from [|12l . |T3 



Tii-i = (1 + exp2Sn-i) _1 is the transmission coefficient of the tunneling process between 
the wells with the least action Sn-i. Depending on the situation of states in each well, 
one of three cases should be chosen in Eq. (|1|): An example of (la) is the decay of high- 
K isomers |J, in which severe breakdown of if-hindrance has been observed recently. 
The present problem of decay-out of near- yrast SD bands corresponds to (lb), where 
the ND states are compound states and described by the GOE model Jlp. In higher 
energies thermally excited SD bands are rotationally damped |7|], and (lc) may be most 
appropriate for the decay-out of such SD continuum states 0. 

In the present case, the relative intensity of 7-rays inside the SD band, which is calcu- 
lated by the total decay-out probability N out (I) at spin /, is almost the only observable 
and will be discussed in the followings. According to [|l[ , iV out is determined by combina- 
tions, r t /D n and r s /r n , of four quantities; the spreading width r t =27iv^/ D n due to the 
coupling ([[]), the level density 1/D n of the ND compound states, and the 7-decay widths 
r s and r n of the SD and ND bands, respectively, where r s is of rotational E2, while r n 
is mainly of statistical El (see 0|| for details). It is worth mentioning that the model 
of IQ for iV out was re-examined by means of a statistical model of compound nucleus 0, 
and both models were found to give identical results for actual range of four parameters 
being relevant to decay of both the A w 150 and 190 SD nuclei. 
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Figure 2. Same as Fig. [I] but for nuclei in the A rs 190 region. Here the back-shift of 1 
MeV is employed throughout in the level density formula. 

In this talk, we present the result of systematic calculations for the relative 7-ray 
intensity as a function of angular momentum in comparison with experimental data. 
We also discuss how the decay-out spin is understood from the calculations. Possible 
directions for improvements of the theoretical framework are suggested further. 



2. RESULTS AND DISCUSSIONS 

Quantities to be calculated are r t , D n , r s , and r n as functions of spin. r s is the usual 
rotational E2 width, while the level density and the statistical El width are given by 

p FG (U) = ( v / ^/48)a-3f/-! eX p2v / aZ7, T n E1 = 0.15 x 2.3 x 10~ n NZA 1/3 (U/a) 5/2 , (2) 

in unit of MeV, where U is excitation energy of the SD band from the ND yrast state, and 



a parameter is taken from an empirical analysis of JTTJ including the shell and temperature 
effects. The decay-out spin J out in the A rs 190 region is generally low so that the back- 
shift is used for A rs 190 nuclei, i.e. U is replaced by U — 1 MeV in pfg, throughout the 
present investigation. As for r%, the least action path in the (£2, 7)-plane is solved and the 
frequency uo s and the action S s - n along it are calculated. Therefore all what we need are 
excitation energies of both SD and ND rotational bands, their quadrupole moments, and 
the potential energy surface and the mass tensor. We have used the Nilsson-Strutinsky 
calculation for potential with the pairing correlations included in the RPA order, and the 
pairing hopping model for mass tensor ||; see || for the detailed formulation. 
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Figure 3. Results with fitting the decay-out spin by scaling the mass tensor (see Table |l|). 
Dotted (dashed) curves denote results with no adjustments (1 MeV back-shift). 

Table 1 



Calculated quantities at the decay-out spin which is fitted by scaling the mass (Fig. [3]). 





Cmass 


lout 


A(eV) 


D n (eV) 


r s (meV) 


r n (meV) 


S s -n(h) 


[/(MeV) 






0.7 


37/2 


390 


630 


1.5 


1.3 


3.3 


2.57 


.52 


152 Dy 


1.6 


28 


2.2 


4.8 


18 


4.7 


5.6 


4.68 


.29 


192 Hg 


2.3 


12 


.11 


69 


.058 


5.2 


6.9 


4.24 


.33 


194 Pb 


2.0 


8 


.061 


134 


.011 


4.3 


7.2 


3.89 


.35 



In Figs. [I] and |^ the results of calculation for relative intensities are shown, where 
those for the lowest band in each parity and signature (four bands) are included in one 
nucleus. The basic characteristics of intensity pattern are reproduced in both the A « 150 
and 190 regions; especially the rapid decrease of transitions at lower spins. Apparently, 
however, the decay-out spin does not agree precisely, and the detailed features, like the 
relative ordering of decay-out for excited bands in one nucleus, or relative difference of 
neighbouring nuclei, are not well described. Comparing Figs. p] and j|, decay-out spins are 
rather well reproduced on average for the A « 150 nuclei, while the calculated Tout's are 
still higher than the observed ones for the A w 190 nuclei, even though the back-shift is 
used for them. Although there are some cases where the yrast SD band decays at higher 
spin than the excited ones in Figs. [I] and |2], it is because the energy ordering is determined 
in the feeding-spin region; the 'real' yrast band survives longest in most cases. 

As is shown in Figs. [I] and 0, the decay-out spin is not precisely reproduced in individual 
cases: What are the reasons? Possible sources of ambiguity in our model may be the level 
density and the mass tensor; we believe that the potential energy surface by the Nilsson- 
Strutinsky method is reliable. As for the change of the level density || it was found that, 
for example, a factor about 10 _1 is enough to fit J out in 152 Dy, while factor 10~ 3 — 10~ 4 
is required for 192 Hg; the latter may be out of the range of allowed ambiguity. The mass 
tensor is calculated by using the Fermi-gas estimate of the number of level-crossings [[J, 
and detailed properties of individual nucleus is not included. Therefore, we try to use a 
mass tensor multiplied by a scaling factor C mass and adjust it so as to fit 7 out . Resultant 
intensities for four selected nuclei are shown in Fig. |3|, and the values of C mass and of 
various quantities at 7 out are summarized in Table [p. Again, the correction factor C mass 
for A fa 150 nucleus is rather small but is considerably larger for A « 190 nucleus. 

In order to see how the decay-out spin is determined by the four relevant quantities, we 
note that N out m J (n/2)(r t / D n ) /(-T s / T n ) holds in realistic situations Two quanti- 
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Figure 5. Comparison of calculated and 
measured number of bands for 143 Eu. 



ties, r t /D n and r s /r n , are plotted as functions of spin in Fig. |j. A out takes appreciable 
value when r t /D n > r s / r n , i.e. two curves crosses in Fig. |, and J out is roughly determined 
by their crossing point. Moreover, the decay-out occurs more rapidly as the crossing angle 
between the two curves gets larger. log 10 (-T s /r , n ) decreases more rapidly in the low-spin 
region because of the behaviour of the rotational E2 transitions of SD bands, and 
therefore the decay-out tends to be more rapid when J out is lower; this behaviour is well 
observed in Figs. [I] and |2]. r s /r n takes similar values for SD bands within the A « 150 or 
190 region, while r t /D n takes variety of values mainly due to the tunneling probability, 
which reflects the different potential landscapes depending on configurations and nuclei. 
This is the main reason why the range of calculated 7 out spreads rather broadly in contra- 
diction to experimentally observed trend. In contrast, observed hindrance factors for the 
decay of ^-isomers take broad range of values in accordance with the results of similar 
calculations of tunneling |]]. The observed trend that J out falls in a rather narrow range 
in both the A »s 150 and 190 regions is not easy to understand. 

Finally, we have recently developed a microscopic model of thermally excited SD states 
that describes both the barrier penetration leading to the decay-out to the ND states, 
and the rotational damping causing fragmentation of rotational E2 transitions || . As an 
example, iVb an d, the number of rotational SD bands in the quasi-continuum 7-ray region 
for 143 Eu is shown in Fig. [j| where solid (dashed) curve is a result using Eq. (lc) (Eq. (lb)). 
Inclusion of the barrier penetration reduces A^and dramatically in the low-spin region and 
the calculation nicely agrees with experimental data |fL4"l , which clearly shows that our 
model is capable for the unified understanding of feeding and decay of the SD bands. 



3. SUMMARY 

Systematic comparison of the calculated and observed relative intensities have been 
performed, and it is found that basic features of the decay-out of SD bands are well 
understood within our theoretical framework. However, the decay-out spin is not precisely 
reproduced in the calculation in individual cases. Especially an average value of 7 out for 
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the A ps 190 nuclei is somewhat overestimated, and the range of calculated J ou t's seems too 
wide in both the A ps 150 and 190 regions. It should, however, be noticed that the average 
value of calculated J out in the A ps 190 region is much smaller than that in the A ps 150 
region by about 10h; the difference in the experiment is even more, about 15H, though. 
As a possible source of improvement a simple scaling of the mass tensor is considered. 
We believe that microscopic calculations of mass tensor based on counting the number of 
level crossings has to be done in order to include detailed nuclear structure effects in each 
SD band. For the detailed comparison, the reproduction of excitation energies of the SD 
bands is very important, which is not enough in the present calculations; for example, that 
of 194 Pb in Table [l] is about 1 MeV higher than the experimentally deduced value. For 
this purpose, the improvement of the potential energy surface may also be necessary, e.g. 
by making use of an another mean-field potential like Woods-Saxon or the one obtained 
by Skyrme Hartree-Fock method. These are remaining future problems. 

As an another possible source of improvement, modification of the compound mixing 
model [T[ may be considered. At higher spins where the excitation energy of SD bands is 
not so high, the ND states to which the SD band couples are not chaotic, and then strong 
reduction of mixing amplitudes may be expected PH| . This reduction leads to a delay of 
decay-out, which hopefully helps to decrease J out especially in the A ps 190. There are, 
however, some ambiguities of the amount of reduction, and a precise modeling of this 
effect is necessary to draw a definite conclusion. It is also an interesting future problem. 

This work is supported in part by the Grant-in-Aid for Scientific Research from the 
Japan Ministry of Education, Science and Culture (No. 12640281). 
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